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Abstract. This paper is a mini-review of the atmospheric 
muon and neutrino flux calculations based upon a recent ver- 
sion of CORT code and up-to-date data on primary cosmic 
rays and hadronic interactions. A comparison of calculations 
with a representative set of atmospheric muon data for mo- 
menta below ^ 1 TeV/c is presented. The overall agreement 
between the calculated muon fluxes and the data provides an 
evidence in favor of the validity of adopted description of 
hadronic interactions and shower development. In particu- 
lar, this supports the low-energy atmospheric neutrino fluxes 
predicted with CORT which are essentially lower than those 
used in current analyses of the sub-GeV and multi-GeV neu- 
trino induced events in underground neutrino detectors. 



1 Introduction 

In recent paper by Fiorentini et al. (2001a) an updated ver- 
sion of FORTRAN code CORT (Cosmic-Origin Radiation 
Transport) has been described and some results of low and 
intermediate energy muon and neutrino flux calculations per- 
formed with CORT have been discussed. A detailed com- 
parison with the muon fluxes and charge ratios measured in 
several modern balloon-borne experiments suggests that the 
atmospheric neutrino flux is essentially lower than one used 
in the current analyses of the sub-GeV and multi-GeV neu- 
trino induced events in underground detectors. Some addi- 
tional results which confirm this conclusion were discussed 
in (Fiorentini et al, 2001b). 

Present work can be considered as an addendum to the 
above-mentioned papers. Its main goal is to provide a com- 
parison of the predicted muon fluxes with a more represen- 
tative set of data from balloon-borne and ground-based ex- 
periments, including rather old (but by no means outdated), 
which cover extensive ranges of muon momenta, ' zenith an- 
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' In this short write-up we hmit ourselves with a selected set of 
muon data at momenta p < 1 TeV/c (some of these data were taken 



gles and atmospheric depths. Such a comparison may be 
useful to obtain further insight into the atmospheric neutrino 
problem. 

2 The CORT code 

Like the earlier versions of CORT (Naumov, 1984; Bugaev 
and Naumov, 1984, 1985, 1986, 1987a,b, 1988, 1989, 1990; 
Naumov, 1993; Bugaev et al, 1998), the new Fortran 90 
code implements a numerical integration of a system of one- 
dimensional kinetic equations describing the production and 
transport of nuclei, nucleons, light mesons, muons, neutri- 
nos, and antineutrinos of low and intermediate energies in 
the atmosphere. It takes into account solar modulation, geo- 
magnetic and meteorological effects, energy loss of charged 
particles, muon polarization and depolarization effects. The 
exact relativistic kinematics is applied for description of par- 
ticle interactions and decays. The new code has a number of 
options and run modes. It is rather flexible and permits fast 
modification of many input data with intrinsic switch keys. 

In order to evaluate geomagnetic effects and to take into 
account the anisotropy of the primary cosmic -ray flux in the 
vicinity of the Earth, CORT uses the method of Ref. (Nau- 
mov, 1984) and detailed maps of the effective vertical cutoff 
rigidities by Dorman et al (1971). The maps are corrected 
for the geomagnetic pole drift and compared with the later 
results reviewed by Smart & Shea (1994) and with the recent 
AMS data the on the proton flux in near Earth orbit (Alcaraz 
et al, 2000a). The interpolation between the reference points 
of the maps is performed by means of two-dimensional local 
B-spline. The Quenby-Wenk relation (Dorman efoL, 1971), 
re-normalized to the vertical cutoffs, is applied for evaluating 

from the recent compilations by Vulpescu et al (1998, 2001) and 
Hebbeker & Timmermans (2001)) and do not consider the under- 
ground and underwater muon data. Additional data can be found in 
the comprehensive reference manual by Grieder (2001), and also in 
Refs. (Bugaev et ai, 1998; Naumov, 1998; Naumov et al, 2000a; 
Sinegovskaya and Sinegovsky, 2001; Naumov, 2001). 
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the effective cutoffs for oblique directions. More sophisti- 
cated effects, Uke the short-period variations of the geomag- 
netic field, Forbush decrease, re-entrant cosmic -ray albedo 
contribution, etc., are neglected. We also neglect the geo- 
magnetic bending of the trajectories of charged secondaries 
and multiple scattering effects. Validity of our treatment of 
propagation of secondary nucleons and nuclei has been veri- 
fied using all available data on secondary proton and neutron 
spectra in the atmosphere (Naumov, 1984; Bugaev and Nau- 
mov, 1984, 1985; Naumov & Sinegovskaya, 2000b, 2001). 

The meteorological effects are included using the Dorman 
model of the atmosphere (Dorman, 1972) which assumes an 
isothermal stratosphere and constant gradient of temperature 
(as a function of depths) below the tropopause. Ionization, 
radiative and photonuclear muon energy losses are treated as 
continuous processes. This approximation is quite tolerable 
for atmospheric depths h < 2 x 10^ g/cm^ at all energies 
of interest (Naumov et al, 1994). Propagation of /x"*" and 
/i^ originating from every source (pion or kaon decay) is 
described by separate kinetic equations for muons with def- 
inite polarization at production. These equations automat- 
ically account for muon depolarization through the energy 
loss (but not through the Coulomb scattering). 



3 Primary cosmic ray spectrum and composition 

In the present calculations, the nuclear component of pri- 
mary cosmic rays is broken up into 5 principal groups: H, 
He, CNO, Ne-S and Fe with average atomic masses A of 
1, 4, 15, 27 and 56, respectively. We do not take into ac- 
count the isotopic composition of the primary nuclei and as- 
sume Z — A/2 for A > 1, since the expected effect on 
the secondary lepton fluxes is estimated to be small with re- 
spect to present-day experimental uncertainties in the abso- 
lute cosmic-ray flux and chemical composition. 

We parametrize the spectra of the H and He groups at 
E < 12Q GeV/nucleon by fitting the data of the bafloon- 
borne experiment BESS obtained by a flight in 1998 (Sanuki 
et al, 2000; Sanuki, 2001). For higher energies (but below 
the knee) we use the data by a series of twelve balloon flights 
of JACEE (Asakimori et al, 1998) and the result of an anal- 
ysis by Wiebel-Sooth et al (1998) based upon a representa- 
tive compilation of world data on primaries. Since the fits 
obtained by Wiebel-Sooth et al (1998) turned out to be very 
close to the combined result of the JACEE experiments, the 
model is caUed "BESSh-JACEE" fit. 

Figure 1 shows a comparison between the BESSh-JACEE 
fit, the data from Sanuki et al (2000); Sanuki (2001) (BESS) 
and Asakimori et al (1998) (JACEE), the fit from Wiebel- 
Sooth et al (1998) (shaded areas) and the data from several 
other experiments (Mason, 1972; Ryan et al, 1972; Seo et 
a/., 1991; Beatty etal, 1993; Ichimura ef aL, 1993; Buckley 
et al, 1994; Diehl et al, 1997; Menn et al, 2000; Aharo- 
nian et al, 1999; Bellotti et al, 1999; Boezio et al, 1999a; 
Alcaraz et al, 2000b, c), performed in different periods of so- 
lar activity. The filled areas in Fig. 1 represent power-type 



parametrizations of the spectra derived by Asakimori et al 
(1998); Ryan et al (1972); Ichimura et al (1993); Aharo- 
nian et al (1999); Alcaraz et al (2000b,c) from the original 
data. The legends indicate the publication dates and (when 
relevant) the years of measurements. 
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Fig. 1. Differential spectra of protons (top panel) and helium nuclei 
(bottom panel) as a function of kinetic energy per nucleon. The 
data points and filled/shaded areas are for the experimental data and 
power-law fits of the data as is shown in the legends. The solid 
curves represent the BESS+JACEE best fit from Fiorentini et al. 
(2001a) (see text for details). 

We assume that the spectra of the remaining three nuclear 
groups are similar to the helium spectrum. This assump- 
tion does not contradict the world data for the CNO and Ne- 
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5 nuclear groups but works a bit worse for the iron group. 
Nevertheless, a more sophisticated model would be unprac- 
tical since the corresponding correction would affect the sec- 
ondary lepton fluxes by a negligible margin. 

In this paper we do not consider the effects of solar mod- 
ulation. Therefore the predicted muon and neutrino fluxes 
are to some extent the maximum ones possible within our 
approach. 

4 Nucleon-nucleus and nucleus-nucleus interactions 

All calculations with the earlier version of CORT were based 
on semiempirical models for inclusive nucleon and meson 
production in collisions of nucleons with nuclei by Kimel' 

6 Mokhov (1974, 1975); Serov & Sychev (1973).2 The 
Kimel'-Mokhov (KM) model is valid for projectile nucleon 
momenta above ~ 4 GeV/c and for the secondary nucleon, 
pion and kaon momenta above 450, 150 and 300 MeV/c, re- 
spectively. Outside these ranges (that is mainly within the 
region of resonance production of pions) the Serov-Sychev 
(SS) model was used. 

Both models are in essence comprehensive parametriza- 
tions of the relevant accelerator data. It is believed that the 
combined "KMh-SS" model provides a safe and model-inde- 
pendent basis to the low-energy atmospheric muon and neu- 
trino calculations. However it is not free of uncertainties. 
For the present calculation, the fitting parameters of the KM 
model for meson and nucleon production off different nu- 
clear targets were updated using accelerator data not avail- 
able for the original analysis (Kimel' & Mokhov, 1974, 1975; 
Kalinovsky et ai, 1985). The values of the parameters were 
extrapolated to the air nuclei (N, O, Ar, C). The overall cor- 
rection is less than 10-15% within the kinematic regions sig- 
nificant to atmospheric cascade calculations. Besides that 
energy-dependent correction factors were introduced into the 
model to tune up the output 7r+/7r~ ratio taking into account 
the relevant new data. 

The processes of meson regeneration and charge exchange 



(7r± + Ai 



+ X etc.) are not of critical importance 



for production of leptons with energies of our interest and 
can be considered in a simplified way. Here we use a proper 
renormalization of the meson interaction lengths, which was 
deduced from the results by Vail et al. (1986) obtained for 
high-energy cascades. 

The next important ingredient of any cascade calculations 
is a model for nucleus-nucleus collisions. Here we consider a 
modest generalization of a simple "Glauber-like" model used 
in (Naumov, 1984; Bugaev and Naumov, 1985). Namely, we 
write the inclusive spectrum of secondary particles c (c = p, 
n, TT^, K"^, K^, . . . ) produced in collisions of nuclei as 



dNAB- 



dx 



AB 



dx 



dx 



+ (1 - Cab) [ZScp + iA- Z)Sc„] S{1 - x). 



^See also Refs. Kalinovsky et al. (1985); Sychev (1999) for the 
most recent versions. 



Here dNf^^^cX I dx is the spectrum of particles c produced 
in a collision of a nucleon (A^ = p, n) with nucleus B, x 
is the Feynman variable, and is the average fraction of 
inelastically interacting nucleons of the projectile nucleus A. 
The term proportional to delta function describes the contri- 
bution of "spectator" nucleons from the projectile nucleus. 

In the standard Glauber-Gribov multiple scattering theory 
the quantity is certainly independent of the type of inclu- 
sive particle c. On the other hand, it depends of the type of 
nucleus collision. Indeed, essentially all nucleons participate 
in the central AB collisions (^^g ~ 1)^ while, according 
to the well-known Bialas-Bleszyhski-Czyz (BBC) relation 



,.incl 



1^ 



for the minimum 



(Bialas et al, 1976), ^^g 
bias collisions. 

To use the above model in a cascade calculation one should 
take into account that nucleons and mesons are effectively 
produced in nuclear collisions of different kind. Namely, the 
contribution from central collisions is almost inessential for 
the nucleon component of the cascade but quite important 
for light meson production. Thus one can expect that effec- 
tively ^^jf > ^^g. We use the BBC relation for nucleon 
production by any nucleus while for meson production we 
P*^*- ^Hc^Air ^ ^' where ^ is a free parameter. Variations of 
this parameter within the experimental limits yield a compar- 
atively small effect to the muon fluxes (except for very high 
altitudes) and inessential (< 6%) effect to the neutrino fluxes 
at sea level (Fiorentini et al., 2001a). Effect of similar vari- 
ations of the parameters ^X'^Vir f^'" other nuclear groups is 
completely negligible. 

5 Numerical results for muons 

Figure 2 displays a comparison of the calculated momen- 
tum spectra of and \i~ for 10 atmospheric depth ranges 
A/ii = 15 - 35, 35 - 65, 65 - 90, 90 - 120, 120 - 150, 150 - 
190, 190 - 250, 250 - 380, 380 - 580, and 580 - 890 g/cm^ 
with the data of two balloon-borne experiments CAPRICE 94 
(Boezio et al, 1999, 2000) and CAPRICE 98 (Hansen et al, 
2001) (the latter experimental data are preliminary). The 
nominal geomagnetic cutoff rigidity is about 0.5 GV (4.5 
GV) and the detection cone is about 20° (14°) around the 
vertical direction with the average incident angle of about 
9° (8°) for CAPRICE 94 (CAPRICE 98). The data points 
in Fig. 2 are the muon intensities at the Flux-weighted Av- 
erage Depths (FAD) defined in (Boezio et al., 1999), while 
the filled areas display the calculated variations of the muon 
spectra inside the ranges A/ij (i = 1, . . . , 10). Namely, 
they are obtained by considering the minimal and maximal 
muon fluxes within each range Ahi. The calculations are 
performed for the conditions of the experiment. The parame- 
ter ^ is taken to be 0.685 (the best value, according to Fioren- 
tini et al (2001a)).'' 



'Here we assume for simplicity that the atomic weight of the 
projectile nucleus is not much larger than that of the target nucleus. 

''it has been shown (Fiorentini et al., 2001a) that the indetermi- 
nation of ^ is only significant for /i < (15 — 20) g/cm^. 
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Fig. 2. Differential momentum spectra of /i^ and for 10 atmospheric depth ranges Ahi (are indicated on the left of each panel). The data 
points are from the CAPRICE 94 experiment (Boezio et a!., 1999, 2000) and from the preliminary analysis of the CAPRICE98 experiment 
(Hansen et al, 2001). The filled areas display the expected variations of the muon fluxes within the ranges Ahi. All the data are scaled with 
the factors indicated at the first and third panels on the right. 




Fig. 3. Differential momentum spectra of positive and negative muons at the ground level. The data points are from the experiments 
CAPRICE 94 (Boezio etai, 1999, 2000; Kremere/a/., 1999) and CAPRICE 97 (Kremer ef a/. , 1999). The data are scaled with a factor 
of 3. The solid (narrow) bands and dotted curves are calculated with CORT as is explained in legends. 
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It is important to note that the thickness of the bands is 
relatively small just for the region of effective muon and 
neutrino production that is in the neighborhood of the broad 
maxima of the muon flux (100 — 300 glcw?). By this is meant 
that, in this region, an error in evaluation of the FAD cannot 
introduce essential uncertainty. Outside the region of effec- 
tive production of leptons, the amplitude of the muon flux 
variations increases with decreasing muon momenta on ac- 
count for the strong dependence of the meson production rate 
upon the depth and the growing role of the muon energy loss 
and decay at ft, > 300 g/cm^. 

Figure 3 shows a comparison of the calculated differen- 
tial momentum spectra of /i+ and /i^ with the most accu- 
rate current data obtained at ground level in the experiments 
CAPRICE 94 (BoezioefflZ., 1999, 2000; Kremeref a/., 1999) 
{h = 1000 g/cm2, = 0.5 GV) and CAPRICE 97 (Kre- 
mer et ai, 1999) (h = 886 g/cm^, = 4.2 GV). The de- 
tection cone in both experiments was the same as described 
above for CAPRICE 94. The calculation are done by using 
the KMh-SS interaction model. Variation of the parameter ^ 
from 0.517 (the BBC value) to 0.710 (an experimental upper 
limit derived from the data on interactions of a particles with 
light nuclei).^ leads to the almost negligible (< 3%) effect in 
the ground-level muon flux. For reference, the spectra calcu- 
lated by using the superposition model for nucleus-nucleus 
interactions are also added. The data below ^ 10 GeV/c 
seem to be precise enough in order to conclude that they dis- 
favor the superposition model. 

As is seen from Figs. 2 and 3, there is a substantial agree- 
ment between the calculations with CORT and the data of 
the CAPRICE experiments within a wide range atmospheric 
depths. This agreement is none the worse than it is for the re- 
cent 3D Monte Carlo calculations (Battistoni et ai, 2001a; 
Engel et ai, 2001; Honda et ai, 2001b; Liu et ai, 2001; 
Poirier et al, 2001; Wentz et al, 2001a,b).^ 

Figure 4 shows expected geomagnetic effect for vertical 
differential momentum spectra of and /i^ and for muon 
charge ratio at sea level. A comparison of these pre- 

dictions with experiment will be discussed elsewhere. 

Figure 5 collects the data on the near-vertical differential 
momentum spectrum of /i+ + ji^ at ground level from many 
experiments performed over the past five decades (Caro et 
al, 1950; Owen & Wilson, 1955; Pine etal, 1959; Pakef a/., 
1961; Holmes et al, 1961; Hayman & Wolfendale, 1962a,b; 
Aurela & Wolfendale, 1967; Baber et al, 1968a; AUkofer 
& Clausen, 1970; Allkofer et al, 1970, 1971; Bateman et 
a/., 1971; Nandi & Sinha, 1972a; Allkofer et al, 1975; Ayre 
et al, 1975; Thompson et al, 1977; Baschiera et al, 1979; 
Green et al, 1979; Rastin, 1984a; Tsuji et al, 1998; Kre- 
mer et al, 1999; Le Coultre, 2001). Only the data for muon 
momenta below ^ 1 TeV/c are shown. 

The compilation, of course, is not exhaustive. In partic- 



See, e.g., Ref. (Kowalski & Bartke, 1981) and references 
therein. 

*See also Ref. (Fiorentini et al, 2001b) for the direct compari- 
son with the FLUKA 3D result. 



ular, it does not involve the low-energy data for large geo- 
magnetic cutoffs, the ground-level data obtained at the al- 
titudes far different from sea level, and also indirect muon 
data from underground measurements.' The spectra calcu- 
lated for /i = 10^ g/cm^, 6^ = 0°, and i?c = with the new 
and old (from Bugaev et al (1998)) versions of CORT are 
plotted. The best-fit spectrum to Nottingham data (Rastin, 
1984a) and the spectrum calculated by Agrawal et al (1996) 
are also shown in bottom panel. 
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Fig. 4. Vertical differential momentum spectra of positive and neg- 
ative muons and muon charge ratios at sea level calculated with 
CORT for several geomagnetic cutoffs Rc. The /i^ spectra on the 
top panel are scaled with a factor of 2. 



^These were discussed in detail by Bugaev et ai (1998). 
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Fig. 5. Near-vertical differential momentum spectrum of /i"*" + /i at ground level (top panel) and its rescaled fragment for p > 10 GeV/c 
(bottom panel). The thickness of the band calculated with CORT reflects the indetermination in the parameter ^. 
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Fig. 6. Muon charge ratio at sea level for near-vertical direction. Geomagnetic cutoffs Rc for different experiments are shown in the legend. 
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In fact, not all experiments listed above measure at sea 
level. Moreover, the experiments were performed in differ- 
ent locations and periods. Therefore the solar modulation 
effects, meteorological and geomagnetic conditions are gen- 
erally different in these experiments. The same is true for the 
detection cones, shielding factors, adopted procedures for the 
momentum spectrum unfolding, and so on. All these factors 
might cause a bias when all the data are compared together 
(Hebbeker & Timmermans, 2001). Here we do not pursue 
any correction procedure (like one suggested by Hebbeker 
& Timmermans (2001)) and plot all the experimental data 
as they are. Instead we show the world average fit derived 
by Hebbeker & Timmermans (2001) from an analysis of (al- 
most) the same data set at p > 10 GeV/c. It can be seen that 
the ground-level experiments are in rather poor agreement to 
one another A large share of the data are arranged outside 
(or even far from) the borderline of the world average fit. 
However the data of the most recent experiments, specifically 
CAPRICE 94 (Ki-emer et ai, 1999) and L3h-C (Le Coultre, 
2001) are placed close to or within these bounds.^ 

Fig. 5 suggests that the muon flux calculated with the new 
CORT is systematically lower than that obtained by Bugaev 
et al. (1998) using the earlier version. The deviation is about 
18% for p = 10 GeV/c and about a few per cent for p > 
80 GeV/c. The main reason of this difference is that the 
BESSh-JACEE primary spectrum below 200 GeV/nucleon 
is essentially lower than that used by Bugaev et al. (1998) 
and in the earlier calculations with CORT. Another reasons 
pertain to the changes in the inclusive cross sections and 
improved description of muon propagation. Below ^ 400 
GeV/c the new spectrum is close to the world average fit 
and to the data from CAPRICE 94 and L3h-C experiments 
as well as to the data from several experiments carried out by 
the BESS Collaboration (Sanuki, 2001; Motoki et al, 2001; 
Sanukief aZ.,2001).'' 

Let us now dwell to the muon charge ratio. It has been 
shown (Fiorentini et al., 2001a) that calculations with CORT 
describe well the modern data on the charge ratio obtained 
at different atmospheric depths. Here we consider only the 
ground-level data for near-vertical directions. The world sur- 
vey is presented in Fig. 6. The data are taken from many 
experiments (Owen & Wilson, 1949; Conversi, 1950; Owen 
& Wilson, 1951; Moroney & Parry, 1954; Filosofo, 1954; 
Fukui, 1955; Holmes etai, 1961; Pak ef aZ., 1961;Coates& 
Nash, 1962; Hayman & Wolfendale, 1962a,b; Kawaguchi et 
al, 1965; Rastin et al., 1965; Baber et al, 1968b; Allkofer et 
al, 1968; Aflkofer & Clausen, 1970; Appleton et a/., 1971; 
Allkofer & Dau, 1972; Nandi & Sinha, 1972b; Burnett et 
al, 1973; Abdel-Monem et al, 1973; Ayre et al, 1973; Bax- 



Note that the L3-FC data on muon flux and charge ratio reported 
in ICRC'27 (Le Coultre, 2001) (see also Refs. (Petersen, 2001; 
Ladron de Guevara, 2001)), are preliminary and may be subject to 
refinement (within the systematic error of 7.7%) when analyses are 
complete. 

'Unfortunately, the tabulated muon data are still unavailable 
from the BESS Collaboration. 



endale et al, 1975; Singhal, 1983; Rastin, 1984b; Stephens 
& Golden, 1987; Grandegger, 1993; Jannakos, 1995; Kriz- 
manic et al, 1995; Boezio et al, 2000; Coutu et al, 2000; 
Brancus et al, 2000; Tsuji et al, 2001; Le Coultre, 2001), 
performed at different geomagnetic locations. According to 
CORT (see bottom panel of Fig. 4), the sea-level charge ra- 
tio is slightly affected by the geomagnetic field. However the 
dispersion of the experimental points is too large and they do 
not show some significant correlation with the geomagnetic 
cutoff. 
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Fig. 7. Differential momentum spectra of muons at sea level for 
several zenith angles and angular bins. The data points are from 
AMH spectrometer (Green etai, 1979). The filled areas display the 
expected variations of the muon fluxes inside the bins AO indicated 
at the left. The solid curves are calculated for average zenith angles 
(6) . All the data are scaled with the factors shown in the parentheses 
at the right. 

Figure 7 shows a comparison with the sea-level data from 
the AMH magnetic spectrometer (Texas A & M and Uni- 
versity of Houston Collaboration) measured within different 
zenith-angle bins A9 (Green et al, 1979). The spectra for 
average zenith angles and expected variations of the spectra 
inside the angular bins are shown by solid lines and filed ar- 
eas, respectively. One can see a good or at least qualitative 
agreement with the data for p < 100 GeV/c and 6 < 80°. 
At higher momenta and at large zenith angles the situation 
is spoiled. However, a comparison between the AMH and 
world survey data for vertical (see Fig. 5) suggests that there 
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Fig. 8. Differential momentum spectra of muons at sea level for 9 — 30° and 75°. The data points are from (Kellogg et al, 1978; Jokisch 
et al, 1979; Tsuji et al, 1998). The filled areas display expected variations of the muon fluxes inside the indicated angular bins. The solid 
curves are for the corresponding average zenith angles. 



is some systematic bias in the AMH experiment above 100 
GeV/c. The abnormal scatter of points in the near-horizontal 
bin is indicative of a flaw in the large-angle measurements. 

Figure 8 displays a comparison between the differential 
momentum spectra of muons at ground level calculated with 
CORT and the data from Brookhaven magnetic spectrom- 
eter (Kellogg et al., 1978), Kiel-DESY muon spectrometer 
located at Hamburg (Jokisch et al, 1979), and Okayama al- 
tazimuthal counter cosmic -ray telescope with a magnet spec- 
trometer (Tsuji et al, 1998). The filled areas display the 
expected variations of the muon fluxes inside the angular 
bins Ae = (26 - 34)° (left panel) and (69 - 81)° (right 
panel). The solid curves are for the corresponding average 
zenith angles of 30° and 75°. The calculations presented 
in Fig. 8 are done without taking into account the geomag- 
netic effects (both the primary spectrum cutoff and muon de- 
flection). These are not important under the conditions of 
Brookhaven and Kiel-DESY experiments but requisite for 
Okayama site.'" 

Generally the geomagnetic effects must decrease the low- 
momentum part of the expected muon flux (see top panel 

'"The Okayama telescope location is (34°40' N, 133° 56' E) and 
vertical geomagnetic cutoff rigidity is about 12.4 GV. 



of Fig. 4). The calculated spectra fit rather well the pre- 
cise Brookhaven and Kiel-DESY data but are in rather bad 
agreement (especially for 6 = 30°) with the Okayama tele- 
scope data of poorer statistics; the geomagnetic corrections 
can only aggravate the disagreement. One can conclude from 
Fig. 8 that there is a systematic flaw in the average zenith an- 
gles measured in the Okayama experiment. 

6 Numerical results for neutrinos 

Let us briefly sketch some results on atmospheric neutrinos 
(AN). Due to geomagnetic effects, the low-energy AN spec- 
tra and angular distributions are quite different for different 
sites of the globe. Figures 9 and 10 display our predictions 
for ten underground neutrino laboratories listed in Table 1 . 

Figure 9 shows the I'e^'Pe, and 17^ energy spectra aver- 
aged over all zenith and azimuth angles. The ratios of the AN 
fluxes averaged over the lower and upper semispheres ("up- 
to-down" ratios) are shown in Fig. 10. As a result of geomag- 
netic effects, both the spectra and the up-to-down ratios for 
the following 6 groups of the labs: SOUDAN + SNO + 1MB, 
HPW, NUSEX + Frejus, Gran Sasso + Baksan, Kamioka, 
and KGF are quite distinct for energies below a few GeV. 
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Table 1. List of some underground laboratories. 



L3b/Det0ctor 


Country 


Geograpitical location 


Geomagnetic location 


Token 


SOUDAN 


USA 


48.00°N 


92.00°W 


58.32° 


331.78° 




SNO 


Canada 


46.80°N 


82.00 W 


57.90 ° 


145.79 ° 




1MB 


USA 


41.72°N 


81.27°W 


52.83° 


346.30° 




HPW 


USA 


40.60°N 


111.00°W 


48.71° 


311.28° 




NUSEX 


Italy 


45.86°N 


6.90°E 


47.24° 


89.09° 




Frejus 


France 


45.1 4° N 


6.69°E 


46.59° 


88.59° 




Gran Sasso 


Italy 


42.45°N 


13.57°E 


42.64° 


94.27° 




Baksan 


Russia 


43.30°N 


42.70° E 


38.06° 


121.64° 




Kamioka 


Japan 


36.42°N 


1 37.31 °E 


26.19° 


204.48° 




KGF 


India 


3.00°N 


78.30°E 


3.25° 


149.36° 





Figure 11 depicts the zenith-angle distributions of Ve, Ve, 
I'f^ and 77^ for Kamioka and Gran Sasso labs calculated with 
CORT using its "standard"(KM+SS) model of particle inter- 
action (see sec. 4) and also the TARGET model foTir/K me- 
son production (including the superposition model for col- 
lisions of nuclei) used by Bartol group (Barr et ai, 1989; 
Agrawal et ai, 1996; Lipari et al, 1998). The TARGET 
model for nucleon production is not included into the CORT 
options yet. Instead the KMh-SS model is used everywhere. 
Below, we refer to CORT switched into the TARGET me- 
son production model as "CORTh- TARGET". For compari- 
son, the result of the recent 3D calculation by Battistoni et 
al. (2000) based on the FLUKA 3D Monte Carlo simulation 
package is also shown. It allows to "highlight" the 3D effects 
which are very dependent on neutrino energy and direction of 
arrival. 
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Fig. 9. Scaled 4t7 averaged fluxes of Ve, Ve, v^i, and //p for ten 
underground laboratories (see table 1 for notation). 
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Fig. 10. Up-to-Down ratios of the Ue,V^, i/^, and 77^ fluxes for ten 
underground laboratories (see table 1 for notation). 



A more quantitative comparison is given by Table 2 which 
tabulates the ratios of thei^e, V^, v^, and fluxes calculated 
with CORTh-TARGET and FLUKA to those with the "stan- 
dard" CORT. The fluxes are averaged over all directions and 
over several energy bins. Table 2 also shows the so-called 
flavor ratio" 



evaluated with CORT, CORTh-TARGET, and FLUKA for the 
same energy bins. 

By using Fig. 1 1 and Table 2, one can conclude that the 
present calculations with CORT for neutrino energies below 
1-2 GeV lead to the AN fluxes which are essentially lower 
than those obtained by Barr et al. (1989); Agrawal et al. 
(1996); Lipari et al. (1998) and by Honda et al. (1990, 1995, 
1996) and those are in current use for many analyses of sub- 
GeV and multi-GeV v induced events in underground de- 
tectors.'^ On the other hand, the new AN fluxes are rather 
close to the those obtained with the earlier versions of CORT 
(Bugaev and Naumov, 1987a,b, 1989, 1990). The main dif- 
ferences (for neutrino energies above 1-2 GeV) are due to the 
new data for the primary cosmic-ray spectrum and composi- 
tion (see Section 3). 



"This quantity roughly represents the ratio of e like to ^ like 
single-ring contained events measured in water Cherenkov detectors 
and the "showers-to-tracks" ratio measured in iron detectors. 

'^A comparison between the earlier AN flux calculations has 
been discussed in detail by Gaisser et al. (1996). 
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Fig. 11. Zenith-angle distributions of Ve, Ve, J^/j, and 77,, for Kamioka and Gran Sasso. Solid curves represent the result of CORT obtained 
with its "standard" interaction model. The dashed curves are the result of CORT obtained with the TARGET model for meson production 
and superposition model for collisions of nuclei. The circles are for the result of calculation by Battistoni et al. (2000) based on FLUKA 3D 
code. The distributions are averaged over the azimuth angle and over the energy bins indicated near the curves. 



Table 2. The ratios of Atv averaged AN fluxes obtained with CORT-l-TARGET and FLUKA to those with CORT and neutrino flavor ratios 
calculated with CORT, CORT-l-TARGET, and FLUKA for Kamioka. 





Average fluxes normalized to CORT 






CORT-l-TARGET 






FLUKA 3D 




(GeV) 
















v^, 


0.1-0.2 


1.64 


1.78 


1.72 


1.73 


1.25 


1.45 


1.41 


1.40 


0.2-0.3 


1.51 


1.68 


1.60 


1.61 


1.23 


1.38 


1.34 


1.33 


0.3-0.5 


1.42 


1.61 


1.48 


1.50 


1.17 


1.28 


1.22 


1.22 


0.5-0.7 


1.34 


1.50 


1.36 


1.38 


1.10 


1.19 


1.11 


1.12 


0.7-1.0 


1.28 


1.41 


1.27 


1.30 


1.05 


1.12 


1.04 


1.05 


1.0-2.0 


1.20 


1.30 


1.17 


1.20 


1.02 


1.05 


0.99 


1.01 


2.0-3.0 


1.10 


1.17 


1.03 


1.08 


0.98 


1.03 


0.96 


0.97 


3.0-5.0 


1.02 


1.07 


0.95 


1.01 


0.99 


1.02 


0.95 


0.97 


5.0-7.0 


0.95 


0.99 


0.91 


0.96 


1.01 


1.05 


0.94 


0.99 


7.0-10.0 


0.91 


0.95 


0.89 


0.93 


1.03 


1.09 


0.95 


1.04 


10.0-20.0 


0.87 


0.91 


0.88 


0.91 


1.07 


1.12 


0.96 


1.02 


20.0-30.0 


0.83 


0.87 


0.87 


0.89 


1.10 


1.15 


0.94 


1.02 





Flavor ratios R,, 




CORT 


CORT-l-TARGET 


FLUKA 3D 


(GeV) 








0.1-0.2 


0.52 


0.51 


0.48 


0.2-0.3 


0.52 


0.50 


0.49 


0.3-0.5 


0.51 


0.50 


0.50 


0.5-0.7 


0.50 


0.50 


0.50 


0.7-1.0 


0.49 


0.50 


0.50 


1.0-2.0 


0.47 


0.49 


0.49 


2.0-3.0 


0.44 


0.47 


0.45 


3.0-5.0 


0.40 


0.43 


0.42 


5.0-7.0 


0.36 


0.37 


0.38 


7.0-10.0 


0.32 


0.32 


0.34 


10.0-20.0 


0.26 


0.26 


0.29 


20.0-30.0 


0.20 


0.19 


0.23 
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The AN angular distributions calculated with CORT and 
with FLUKA 3D code (Battistoni et al, 2000, 2001b) are in 
good agreement for any zenith angle at energies above 0.7- 
0.8 GeV. On the other hand, CORT is systematically lower 
than FLUKA for < 1 GeV. At the energy = 0.5 GeV 
the FLUKA An averaged i/^, v^, and fluxes exceed 
the corresponding CORT results by about 15%, 27%, 18% 
and 18%, respectively. This discrepancy is only partially due 
to 3D effects (Fig. 11), which can account for an increase 
< 10%, and it is most probably related to differences be- 
tween the hadronic interaction models adopted in FLUKA 
and CORT codes. 

It is pertinent to note that several recent 3D Monte Carlo 
calculations of the AN flux (Honda et al, 2001a,b; Liu et al., 
2001; Ply askin, 2001; Tserkovnyakef a/., 2001; Wentz et al, 
2001a) confirmed the sharp enhancement of low-energy neu- 
trino intensities for near-horizontal directions predicted with 
FLUKA. However the quantitative characteristics of this en- 
hancement are very model-dependent and vary over a broad 
range from one calculation to another 

Prompt neutrinos 

Figures 12 and 13 collect the differential energy spectra of 
downward going atmospheric neutrinos calculated within a 
wide energy range (from 50 MeV to about 20 EeV) for 1 1 
zenith angles. Figures show the "conventional" neutrino con- 
tribution (originated from decay of pions, kaons and muons) 
and the total AN spectra which include the "prompt" neutrino 
contribution originated from semileptonic decays of charmed 
hadrons (mainly D^, D^, D mesons and hyperons).'"* 
The prompt neutrino contribution must dominate at very high 
energies. However the charm hadroproduction cross sections 
are very model-dependent and cannot be unambiguously pre- 
dicted for lack of a generally accepted model. As a result 
the prompt neutrino contribution and even the energies above 
which the prompt muon and electron neutrinos become dom- 
inant are very uncertain as yet (Bugaev et al, 1998; Pasquali 
et al, 1998, 1999; Gelmini et al, 2000a,b; Volkova & Zat- 
sepin, 2001; Costa, 2001; Costa et al, 2001). The results 
shown in Figs. 12 and 13 are obtained by using the two 
phenomenological approaches to the charm production prob- 
lem: the quark-gluon string model (QGSM) and recombina- 
tion quark-parton model (RQPM). The prompt muon fluxes 
predicted by QGSM and RQPM are both consistent with the 
current deep underground data and may be considered as the 
safe lower and upper limits for the prompt muon contribu- 
tions. The basic assumptions of these models and the as- 
pects of atmospheric cascade calculations were described by 
Bugaev ef a/. (1989). 

'^We do not discuss here the contributions from Da mesons and 
heavier states (such as bb, tt, W, Z) which become important 
sources of atmospheric (and extraterrestrial) prompt r neutrinos at 
very high energies (Pasquali & Reno, 1999; Athar et al, 2001). 

'''See also Refs. (Bugaev et al., 1998; Naumov, 1998) for more 
details and related references. 



Calculations of conventional AN fluxes for energies below 
^ 100 GeV are done with CORT for Kamioka site, while the 
high-energy part of the conventional neutrino spectra as well 
as the prompt neutrino contributions are obtained by using 
the results of Ref. (Naumov et al, 1998). Needless to say 
the high-energy AN fluxes are independent of location. The 
high-energy calculation takes into account many "thin" ef- 
fects, like A'e3 and K^-^ form factors and K'g semileptonic 
decays, meson regeneration and charge exchange through re- 
actions 7r± + Air ^ 7r±(=F) + 7r± + Air ^ K^^°^ + X, 
etc., as well as through Ktj2 and K^,^ decays. The low- and 
high-energy parts of the spectra are smoothly merged by us- 
ing a polynomial (over \og{E^)) clutching functions. The 
results of these calculations for all underground detectors are 
now embedded in a simple FORTRAN code "CORTout". It 
is based on a two-dimensional spline interpolation over the 
detailed tables and it may be useful for a fast evaluating the 
energy spectra and zenith-angle distributions of conventional 
and prompt atmospheric neutrinos within the energy range 
from 0.05 to ^ 10^° GeV. The code is available upon request 
from the author 

7 Conclusions 

Let us briefly re-state the main conclusions which follow 
from the outcome of Fiorentini et al. (2001a,b) and from the 
present study. 

The results of extensive calculations performed with the 
new version of code CORT by using the updated KMh-SS 
hadronic interaction model and the BESSh-JACEE fit for the 
spectrum and composition of primary cosmic rays are in sub- 
stantial agreement with a quite representative set of data on 
muon momentum spectra and charge ratios measured at dif- 
ferent atmospheric depths, zenith angles, and geomagnetic 
locations. This provides a solid evidence for the validity of 
our description of hadronic interactions and shower devel- 
opment. The low-energy atmospheric neutrino fluxes calcu- 
lated with CORT for many geomagnetic locations are sys- 
tematically lower than those used in the current analyses of 
the data from underground neutrino experiments, while the 
neutrino flavor ratio and "up-to-down" asymmetry are essen- 
tially the same as in all recent calculations. It is difficult to 
increase the AN flux without spoiling the agreement with the 
current data on hadronic interactions, primary spectrum and 
muon fluxes. 
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Fig. 12. Differential energy spectra of electron neutrinos and antineutrinos for 11 zenith angles 9. Low-energy range is for Kamioka site. 
Contributions from charm decay are calculated using recombination quark-parton model (RQPM) and quark- gluon string model (QGSM). 
At high energies, from smallest to largest fluxes, cos 6 varies from to 1 with an increment of 0. 1 for each group of curves. 
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Fig. 13. Differential energy spectra of muon neutrinos and antineutrinos for 11 zenith angles 9. Low-energy range is for Kamioka site. 
Contributions from charm decay are calculated using recombination quark-parton model (RQPM) and quark-gluon string model (QGSM). 
At high energies, from smallest to largest fluxes, cos 6 varies from to 1 with an increment of 0. 1 for each group of curves. 
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